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l> ■ 

■ Evolution of anisotropy in momentum and coordinate space of the parton system produced in 
\ relativistic heavy-ion colhsions is discussed within the free-streaming approximation. The momen- 

^vq . tum distribution evolves from the prolate shape - elongated along the beam - to the oblate one - 

^ ■ squeezed along the beam. At the same time the eccentricity in the coordinate space, which occurs 

finite values of impact parameter, decreases. It is argued that the parton system reaches local 
, thermodynamic equilibrium before the momentum distribution becomes oblate. 

l/^ ■ PACS numbers: 25.75.-q,12.38.Mh 

I. INTRODUCTION 

^ , Parton system, which emerges at the early stage of relativistic heavy- ion collisions, is anisotropic both in momentum 
^ ■ and coordinate space. The anisotropics crucially influence the dynamics of the system: the momentum one causes 
^ I plasma color instabilities (for a review see [1]); the coordinate-space one is responsible for the hydrodynamic elliptic 
flow (for a review see [2]). An eccentricity of the overlap region of colliding nuclei at nonzero impact parameter 
decreases when the parton system produced in the overlap region expands. While the eccentricity simply decays, the 
parton momentum distribution, which is observed locally, changes from a strongly prolate shape - elongated along 

■ the beam axis - to an oblate form - squeezed along the beam. 

\ Since hydrodynamics requires at least partial local equilibrium Q, an observation of the elliptic flow suggests that 
CS| ■ the parton system is thermalized before the initial eccentricity is significantly reduced. The equilibration time of the 

\ parton system teq was actually estimated to be shorter than 1 fm/c [4]. Inter-parton collisions cannot equilibrate the 
\^ . system so fast but magnetic unstable modes due to the momentum anisotropy speed-up the equilibration process. 

' However, the question arises whether the momentum distribution is prolate or oblate just before the equilibrium is 

. reached. In this note we attempt to resolve the issue in a very simple classical model where partons produced in the 

■ overlapping region of colliding nuclei freely escape from it. We analyze how the coordinate space anisotropy decays 
and how the momentum distribution evolves in a box which includes the Lorentz contracted region where the partons 
are initially produced. The effect of finite formation time of produced partons is taken into account. 

We are fully aware how naive our approach is. Matter, which emerges at the early stage of relativistic heavy- ion 
^ . collisions, is very dense, presumably inhomogeneous, and partially coherent due to the memory of pure quantum 
state of two colliding nuclei. Such a system cannot be reliably described in terms of kinetic theory with weakly 
interacting quasi-particles on mass-shell. A derivation of transport equation from quantum field theory clearly reveals 
the limitations of the kinetic approach (sf. We believe, however, that our free-streaming model still grasps global 
features of the early stage system. The evolution of anisotropics, which is of our main interest here, is dominated 
by the system's expansion. And it proceeds with the velocity of light independently of details of the system's 
dynamics. However, it should be clearly understood that in the free-streaming model, where no interaction is present, 
a real equilibration does not take place. Therefore, even if the particle momentum distribution appears to be of 
the equilibrium form, one cannot conclude that the system has reached the equilibrium, which is required by the 
hydrodynamic description, as the pressure resulting from inter-particle collisions is absent in the free-streaming model. 
We return to the discussion in the concluding section of the paper. 



X 



* Electronic addres: weronika. jasOgmail . com 
^ Electronic address: mrowOf uw . edu . pi 



2 



II. DECAY OF ECCENTRICITY 



The partons are assumed to be produced in an ellipsoidal region parameterized by the three-dimensional Gaussian 
function centered at zero with the widths: (Jy^ ctz^ where y and z denote Cartesian coordinates. As usually, the 
axis z is along the beam. The distribution of parton's rapidity, which is denoted here by F, is assumed to be Gaussian 
as well with the width AK. We note that the rapidity distribution of charged pions produced in Au-Au collisions 
at y/s = 200 GeV per nucleon-nucleon pair is well described by the Gaussian distribution with AK = 2.3 [6]. Since 
we work in the center-of-mass frame, the rapidity distribution is centered at zero. We also assume that partons are 
massless and then, as it will be evident later on, we do not need to specify the distribution of their transverse momenta 
which is denoted as P{pt)- Thus, the distribution function of partons, which obeys the collisionless Boltzmann 
equation, is 
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with (j) being the azimuthal angle in the momentum space. We note that we use the natural units where c = 1. The 
distribution function is normalized to unity 
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When the impact parameter of the colliding nuclei is chosen to be along the axis as shown in Fig. [TJ the 
eccentricity, which drives the elliptic flow, is defined as 
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With the distribution function ([T]), one computes 
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If one takes into account only partons with the rapidities Y obeying Fmin <Y < Fmax, the formula, which defines the 
coefficient a, changes to 

jy""" dYchr'^Y exp [ - 1^2] 
a = — 2iiiL_ _ ^ 
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Computing {y'^) analogously to (x^), one finds the eccentricity ([2]) as 
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where the initial eccentricity €{0) and the average transverse size of the overlap region of colliding nuclei Rt are 
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FIG. 1: (Color online) A view of colliding nuclei of equal 
radius as seen in the plane x—y transverse to the beam 
axis. 



FIG. 2: (Color online) The eccentricity in Au-Au collisions 
as a function of time for five values of the impact parameter. 
The curve with the largest initial eccentricity corresponds to 
b = 1.8R, the next one to 6 = 1.6R, etc. 



The formula ([3]) was earlier derived in [H, 0] for a narrow interval around Y = when a = 1/2. Unfortunately, by 
mistake a = 1 in [7]. 

When the impact parameter varies, both the initial eccentricity e{0) and the average transverse size Rt change. 
To express the two quantities through the nuclear radius R (a hard sphere parameterization is adopted here) and the 
impact parameter 6, we replace the overlap region of two circles shown in Fig. [T] by the ellipse given by the equation 

with the half- axes dx and dy defined in Fig. [TJ Elementary geometric arguments provide dx = R — b/2 and dy = 
^JR? — 1^. Since the mean square sizes of the ellipse are (x^) = and (^^) = we identify and Oy with 

dxi'l and dyj2^ respectively. Then, one expresses £:(0) and Rt through R and h as 

In Fig. [2] we show predictions of the formula (j3|) for Au-Au colhsions at midrapidity (a = 1/2). The radius of a 
gold nucleus is chosen to be 7 fm. As seen in Fig. O the larger initial eccentricity, the faster its decay. We note that 
the largest elliptic flow in Au-Au collisions is observed at 6 ~ 10 fm Q corresponding to 6 ~ \AR. At larger impact 
parameters, the produced system is presumably too small to fully manifests collective hydrodynamic behavior. An 
analysis of the experimental elliptic flow data within the hydrodynamic model shows that the eccentricity cannot be 
reduced to more than 75% of its initial value [1], see also 0. When the reduction is larger, the ideal hydrodynamics, 
which gives an upper limit of the flow, significantly underestimates experimental data. Fig. [2] shows that for h ^ lAR 
the eccentricity is reduced to 75% of its initial value at 1.5 fm/c. Within this time interval the system has to be 
equilibrated to start hydrodynamic evolution responsible for the elliptic flow. Thus, our estimate of the upper limit 
of equilibration time is 1.5 fm/c. Actually, the hydrodynamic analysis [4], which uses not only the elliptic flow data 
but other experimental constraints as well, provides the equilibration time teq as short as 0.6 fm/c. 



III. MOMENTUM ANISOTROPY EVOLUTION 



In this section we compute the momentum distribution of partons in a box of sizes L^^ Ly and L^. The box is 
centered dit x = y = z = 0. To simplify the calculations, the sharp edge box is replaced by the 'box' function of 
Gaussian form 
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FIG. 3: (Color online) The momentum anisotropy as a 
function of time for four values of the rapidity distribution 
width /Sy. The most upper line corresponds to = 2.5, 
the lower one to AK = 2.0, etc. The longitudinal size of 
the box Lz equals 1 fm. 



FIG. 4: (Color online) The momentum anisotropy as a 
function of time for four values of the rapidity distribution 
width /Sy. The most upper line corresponds to = 2.5, 
the lower one to = 2.0, etc. The longitudinal size of 
the box Lz equals 3 fm. 



which obeys the conditions 
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The momentum distribution of particles in the box is characterized by the parameter 
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with analogous formula for /(t,r, p) is given by Eq. ([T]). For the prolate distribution, one has p > 1, for the 

isotropic one, there is p = 1, and finally, for the oblate distribution, we have p < 1. Since the parameter p is sensitive 
only to the second moments of the momentum distribution, p = 1 is a necessary but not a sufficient condition for an 
isotropy of the distribution. 

The calculations simplify when the system is cylindrically symmetric in coordinate space that is (Jx = (Jy = cft and 
Lx = Ly = Lt- Then, the anisotropy parameter ([5j) is given by 
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One easily computes the initial value of p as p{Q) = e^'^^ — 1. It is also of interest to see how the energy density e{t) 
in the box decreases when the momentum anisotropy evolves. A simple calculation provides 
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To compute p{t) and e(t)/e(0), the parameters cr^, cr^? Lt and have to be chosen. As well known, the nuclear 
density of heavy nuclei is well described by the Woods-Saxon formula which can be roughly approximated by the 
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FIG. 5: (Color online) The relative energy density as a 
function of time for four values of the rapidity distribution 
width /Sy. The most upper line corresponds to = 1.0, 
the lower one to AK = 1.5, etc. The longitudinal size of 
the box Lz equals 1 fm. 



FIG. 6: (Color online) The relative energy density as a 
function of time for four values of the rapidity distribution 
width /Sy. The most upper line corresponds to = 1.0, 
the lower one to = 1.5, etc. The longitudinal size of 
the box Lz equals 3 fm. 



sharp-sphere parameterization with the radius R which for the heaviest nuclei equals about 7 fm. The mean square 
radius for the sharp-sphere parameterization equals (r^) = 3R^/5. Therefore, we choose the widths of the Gaussian 
distribution (Tt = (Tx = ay to be equal to R/y/b ~ 3fm. The transverse size of the box Lt is assumed to coincide 
with (Jt- The longitudinal width of the interaction zone cjz is chosen as 1 fm. The calculations of p{t) and e(t)/e(0) 
are performed for = 1 fm and 1/^ = 3 fm. The results are shown in Figs. O [H [5]and[6l 

As already mentioned, the width of rapidity distribution of charged pions produced in Au-Au collisions at ^/s = 200 
GeV is /Sy = 2.3 [61]. One expects that the rapidity of produced partons is even broader. Therefore, the results shown 
in Figs. [3l [Hfor AK = 2.5 seem to be relevant for heavy-ion collisions at RHIC. In such a case, it takes 6-8 fm/c to 
have an oblate momentum distribution. As Figs. [5l [6] show, the energy density in the box is then decreased by a large 
factor - the system is much diluted. 

One argues that particles, which are produced with rapidity F, materialize only at a finite proper time r and 
space-time rapidity rj = Y. Keeping in mind that 



r = \J — , T] = -In 

one finds 



t = r chf] , z = r shr] . 

When the formation time r is 0.3 fm/c and r] = Y = 2.5, one obtains t = 1.8 fm/c and z = 1.8 fm. Thus, partons 
with y = 2.5 materialize beyond the box of L;^ < 1 fm. To take into account the effect of finite time formation in our 
analysis, we simply eliminate from the distribution function ([T]), the partons of the proper time smaller than r. In 
other words, the function ([T]) is multiplied by Q{t^/l — v1 — r) and high-rapidity partons are effectively excluded. 

The temporal evolution of the momentum asymmetry p, which takes into account the finite formation time, is 
shown in Figs. [3 and [8] for r = 0.3 fm/c and r = 0.8 fm/c, respectively. As previously, = Lt = S fm, = 1 fm 
and Lz = 3 fm. As seen, the effect of finite formation time is very significant. At the beginning the momentum 
distribution is oblate as partons with small Y appear in the box earlier than those with larger Y. After some time 
the momentum distribution is prolate, and it becomes again oblate due to the system's expansion only after 3-5 fm/c. 



IV. CONCLUSIONS AND DISCUSSION 



As discussed in Sec. [Ill the parton system produced in nucleus-nucleus collisions has to be equilibrated before 
1.5 fm/c. Otherwise the eccentricity is too much reduced and the ideal hydrodynamics significantly underestimates 
experimental data of the elliptic flow. The results from Sec. Illll show that at the time 1.5 fm/c the local momentum 
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FIG. 7: (Color online) The momentum anisotropy as a 
function of time for the formation time r = 0.3 fm/c and 
four values of the rapidity distribution width AK. The 
most upper line corresponds to /SV = 2.5, the lower one to 
= 2.0, etc. The longitudinal size of the box Lz equals 
3 fm. 



FIG. 8: (Color online) The momentum anisotropy as a 
function of time for the formation time r — 0.8 fm/c and 
four values of the rapidity distribution width /SV. The 
most upper line corresponds to /SV = 2.5, the lower one to 
= 2.0, etc. The longitudinal size of the box Lz equals 
3 fm. 



distribution is still prolate. Therefore, we conclude that just before local equilibrium is reached, the parton momentum 
distribution is elongated along the beam. Let us consider how reliable is the conclusion. 

First of all, if our rather conservative estimate of the upper limit of equilibration time of teq = 1-5 fm/c is replaced 
by more elaborated estimate of 0.6 fm/c given in [1], our conclusion is really safe - it seems impossible to build up 
the oblate momentum distribution in such a short time. 

One wonders whether the formation time r can be extended. It should remembered, however, that in the case of 
finite r, the initial coordinate-space eccentricity does not occur at t = but at t = r (for mid rapidity particles). And 
the interval of time when the system reaches equilibrium is reduced to tgq — r. So, the longer r, the fast thermalization 
is more and more difficult to understand. 

Temporal evolution of the momentum anisotropy is faster when a^, which is the initial longitudinal localization 
of produced partons, is reduced. However, for = 3 fm, the time when p = 1 is shorter only by 25% when 
decreases from 1 fm to 0.5 fm. So, our conclusion remains unchanged. We note that should not be confused 
with the longitudinal localization of valence quarks of incoming nuclei which, due to the Lorentz contraction, is (Jt/t 
with 7 being the Lorentz factor. Our corresponds to the produced partons. Therefore, it cannot be too small, 
as the wee partons, which are localized beyond the contracted volume of incoming nuclei, effectively participate in 
nucleus- nucleus collisions. 

Our free-streaming model of massless partons is obviously very naive. Partons, which are produced at the collision 
early stage, often carry a large virtuality acting as a mass. If the parton mass is taken into account both the 
coordinate and momentum space evolutions are slowed down. Inter-parton interactions presumably lead to a similar 
effect. However, we cannot see a good reason that the coordinate space evolution is slowed down much more than the 
momentum space evolution. Therefore, our conclusion seems to be rather safe. 

Obviously it is desirable to improve our free-streaming model but the problem is rather difficult. Color Glass 
Condensate (CGC) approach (for a review see [10]), which is the best developed effective theory to study the early 
stage of heavy-ion collision, is not well suited for the problem. The valence quarks are treated in CGC as classical 
passive sources of small x, highly populated gluons which are described in terms of classical fields. Partons with 
sizeable x, which crucially influence the momentum distribution, are essentially absent in CGC. 

Numerous studies of equilibration of parton systems with the initially oblate momentum distribution, which are 
reviewed in [Ij], are theoretically well founded as the oblate system is diluted and decohered. However, according to 
our conclusion, these studies do not actually explain how local equilibrium is reached in heavy-ion collisions but rather 
how the equilibrium is sustained. We believe that the plasma equilibration, in particular the color instabilities, which 
are supposed to speed up the process of thermalization, should be analyzed as in [ll|: [HI ^^^^ systems with 

prolate momentum distribution. We are aware that kinetic theory is hardly applicable to such dense, inhomogeneous 
and partially coherent systems but Nature does not seem to be bothered by our theoretical difficulties. 



7 



Acknowledgments 

We are grateful to Uli Heinz for helpful correspondence. This work was partially supported by Polish Ministry of 
Science and Higher Education under grant 1 P03B 127 30. 



[1] St. Mrowczynski, Acta Phys. Polon. B 37, 427 (2006). 

[2] P. F. Kolb and U. W. Heinz, in Quark Gluon Plasma 3, edited by R.C. Hwa and X.N. Wang (World Scientific, Singapore, 
2004). 

[3] P. Arnold, J. Lenaghan, G. D. Moore and L. G. Yaffe, Phys. Rev. Lett. 94, 072302 (2005). 

[4] U. W. Heinz, AIP Conf. Proc. 739, 163 (2005). 

[5] St. Mrowczynski and P. Danielewicz, Nucl. Phys. B 342, 345 (1990). 

[6] I. G. Bearden et al [BRAHMS Collaboration], Phys. Rev. Lett. 94, 162301 (2005). 

[7] P. F. Kolb, J. Sollfrank and U. W. Heinz, Phys. Rev. C 62, 054909 (2000). 

[8] C. Adler et al [STAR Collaboration], Phys. Rev. C 66, 034904 (2002). 

[9] H. J. Drescher, A. Dumitru, C. Gombeaud and J. Y. Ollitrault, arXiv:0704.3553 [nucl-th]. 

[10] E. lancu and R. Venugopalan, in Quark Gluon Plasma 3, edited by R.C. Hwa and X.N. Wang (World Scientific, Singapore, 
2004). 

[11] St. Mrowczynski, Phys. Rev. C 49, 2191 (1994). 

[12] J. Randrup and St. Mrowczynski, Phys. Rev. C 68, 034909 (2003). 



